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Abstract 

We study Lagrangian time-discretizations of the Hunter-Saxton equa- 
tion. Using the Moser-Veselov approach, we obtain such discretizations 
defined on the Virasoro group and on the group of orientation-preserving 
diffeomorphisms of the circle. We conjecture that one of these discretiza- 
tions is integrable. 

2000 Mathematical Subject Classification 34K99, 22E65, 70H99 
Keywords: Hunter-Saxton equation, Virasoro group, group of 
orientation-preserving diffeomorphisms of the circle, discrete Lagrangian 
system 



Introduction 

In 1991 J. K. Hunter and R. Saxton considered the equation 

(itt + UU:c)x = ^ul, (1) 

which describes the propagation of weakly nonlinear orientation waves in 
a massive nematic liquid crystal director field. A derivative with respect 
to X of the Hunter-Saxton equation 

1 2 

{ut + uux)xx = {-y-x)x, or simply utxx = —2uxUxx — uuxxx, (2) 

is also often called the Hunter-Saxton equation. It is equation (^) that 
will be central to the present paper. 

Let us recall some results concerning the Hunter-Saxton equation. 
Equation (jl|) was solved in ^ using the method of characteristics. A 
generalization of equation (hi) was studied by M. V. Pavlov || and it was 
also solved. Equation was investigated by J. Hunter and Y. Zheng [^, 
and it was proven that (^) is a completely integrable, bi-variational, bi- 
Hamiltonian system. 

The Hunter-Saxton equation is related to the Korteweg-de Vries and 
Camassa-Holm equations in many ways. P. Olver and P. Rosenau de- 
scribed Q this relation using tri-Hamiltonian structures and scaling ar- 
guments. R. Beals, D. H. Sattinger and J. Szmigielski described ^ the 
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scattering theory for all three equations in a unified way. B. Khesin and 
G. Misiolek interpreted these three equations as the Euler equations 
describing geodesic flows associated to different right-invariant metrics on 
the Virasoro group, or on an appropriate homogeneous space. Moreover, 
they proved that these three equations exhaust all generic bi-Hamiltonian 
systems related to the Virasoro group. The paper was a motivation 
for the present paper and we will recall some results of in Section ^. 

Since the Hunter-Saxton equation ^ is an Euler equation, it is natural 
to try to discretize it using the Moser-Veselov approach. This approach 
is the following. 

Let Ad be a manifold, let L be a function on Ad x Af. A discrete 
Lagrangian system (see the review |^ for references) describes stationary 
points of a functional S = S{X) defined on the space of sequences X = 
{xk),Xk G AI, fc G Z, by a formal sum 

S{X) = ^ L{xk,Xk+i). 

feez 

The function L is called the Lagrangian. 

Let us assume that M is a Lie group G. It was observed by A. P. Veselov 
and J. Moser that symmetric 

L{x,y) = L{y,x) 

and right-invariant (or left-invariant) 

VffGG L{xg,yg) = L{x,y) 

Lagrangians often correspond to integrable systems which are discretiza- 
tions of the Euler equations corresponding to some right-invariant (or 
left-invariant) metrics on G. The case of SO{n) and other classical finite- 
dimensional groups U{n), Sp{n) was studied extensively in ^, 
The first example of an infinite-dimensional group, namely SDiff(R ) was 
considered in |l^ . 

We shall study discrete Lagrangian systems on the Virasoro group 
in order to find discretizations of the Hunter-Saxton equation. Dis- 
crete Lagrangian systems on the Virasoro group were studied in jl^, |l^ . 
These papers were motivated by the observation by V. Yu. Ovsienko and 
B. A. Khesin jl6j that the Korteweg-de Vries equation is the Euler equa- 
tion for a right-invariant metric on the Virasoro group. In [l4| [l^ we 
investigated discrete Lagrangian systems on the Virasoro group trying to 
find a time-discretization of the Korteweg-de Vries equation. In we 
studied a particular class of discrete Lagrangian systems containing an in- 
teresting system with remarkable properties. It was noticed however that 
this system had no visible relation to the Korteweg-de Vries equation. We 
will show in the present paper that this system is in fact related to the 
Hunter-Saxton equation. In |l5j we considered another particular class of 
discrete Lagrangian systems which has the Korteweg-de Vries equation as 
a continuous limit. One of these systems has very interesting properties 
and we conjectured that this system is integrable, but in fact till now only 
one integral was found. 
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In this paper we study a class of discrete Lagrangian systems on the 
Virasoro group satisfying natural conditions generalizing those considered 
in 1^^. We show that these systems have the Hunter-Saxton equation as 
their continuous limit. We show that the system considered in which 
belongs to this class can be considered as a good candidate for a natural 
integrable Lagrangian time-discretization of the Hunter-Saxton equation. 

It is shown in that the Hunter-Saxton equation can also be consid- 
ered as an Euler equation on the homogeneous space Diff+(S'^)/Rot(S'^) 
of the group of orientation-preserving diffeomorphisms of the circle mod- 
ulo the group of rotations of the circle. We will also consider possible 
discretizations of the Hunter-Saxton equation on the group Diff+(5'^). 
However, our primary interest will be the Virasoro group since it plays 
a more fundamental role for the Hunter-Saxton equation. Moreover, we 
will see that in some sense discrete Lagrangian systems on Diff+(5'^) are 
"reductions" of discrete Lagrangian systems on the Virasoro group. 

The plan of the paper is as follows. In Section ^ we recall necessary 
facts about the Virasoro group and the Virasoro algebra. In Section H 
we recall the results of paper Q describing links between the Hunter- 
Saxton equation and the Virasoro group. In Section plwe study discrete 
Lagrangian systems on the Virasoro group. In Sectionkl we recall the re- 
sults of |l4[ and show that the system studied in can be considered as 
a good candidate for a natural integrable Lagrangian time-discretization 
of the Hunter-Saxton equation. In the final Section ^ we will study dis- 
crete Lagrangian systems on the group Diff+(S'^) and their relation to the 
Hunter-Saxton equation. We will see that the "reduction" to Diff+(S'^) 
of the system studied in Section ^ is particularly simple. 

1 The Virasoro group and the Virasoro 
algebra 

Let Diff+(S^) be the group of diffeomorphisms of S preserving the ori- 
entation. We shall represent an element of Diff+(S'^) as a diffeomorphism 
/ : R ^ K such that 



Such a representation is not unique. Indeed, the functions / + 2nk, k £ Z 
represent one element of Diff+(5'^). 

There exists a non-trivial central extension of Diff+(S'^) which is 
unique up to an isomorphism 

This extension is called the Virasoro group (or the Bott- Virasoro 
group) and is denoted by Vir. Elements of Vir are pairs {f,F), where 
/ £ Diff+ {S^ ) , F £ R. The product of two elements is defined with the 
help of the Bott cocycle as 



3, 



2 



1 



/ G C°°(R), 
f'(x) > 0, 

fix + 27t) = fix) + 27V. 



if,F)o{g,G) = {fog,F + G+ / log(/ o g)' d log <?')• 
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The unit element of Vir is (id, 0). The inverse element of {f,F) is 
{f-\-F). 

The Virasoro algebra vir is a Lie algebra corresponding to the Virasoro 
group. It is the central extension of the algebra vect{S^) of vector fields 
on the circle 5^ 

vir = vect{S^) R. 

We represent an element of vect{S^) as v{x)dx, where w is a 27r-periodic 
function. Thus an element of the Virasoro algebra is a pair {v{x)dx,a)- 
The algebra commutator in vir is defined with the help of the Gelfand- 
Fuchs cocycle as 

[{v{x)dx,a], [w{x)dx,b]) = {{—vwx + Vxw){x)dx, J v^xxW dx). 





2 The Hunter-Scixton equation and the 
Virasoro group 

In this section we briefly recall some results of |Q| concerning a link between 
the Hunter-Saxton equation and the Virasoro group. 

Let G be a Lie group, g its Lie algebra. A right-invariant metric on G 
is completely defined by its restriction to g 

E{v) = ^{v,Av), uGg, 

where ^4 is a linear map A : g — > g* called the inertia operator. 

To describe a geodesic g{t) on G, we transport its velocity vector to 
the identity by the right translation 

Since v{t) is an element of g, we can consider m — Av which is an element 
of the dual space g*. Then m satisfies the Euler equation given by the 
following explicit formula: 

dm , 

— = -adj^-i^m. (3) 

This is a standard result which can be found in If we start with a 

left-invariant metric, the sign in is reversed. 

Let us now compute ad* for the Virasoro algebra. Let a space 

{{u{dx)^ , a)\u G ,u(x) = u(x + 2-k), a G R} 

be the dual space vir* to the Virasoro algebra with the natural pairing 
given by 

{{u{dx)'^ ,a), (wdx, c)) = J uw dx + ac. 
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Thus, ad* is defined by the formula 

(arf(^s,,6)("('^a;)^a), (wa:r,c)) = {{u{dxf ,a),[{vd^,b), [wd^, c)]) . 

Using the definition of the commutator in vir and integration by parts we 
obtain that the right-hand-side is equal to 





Thus we obtain 



Z7r 



a) = {{2uv^ + u^v - av^.j,x){dxf , 0). 

Now let us consider an inertia operator A : vir — > vir*, given by 

A(«a.,6) = ((-A^;)(dx)^&), 

where A = —d%. The corresponding scalar product on g is given by the 
formula 



E((vdx,b)) = — y vvxx dx + 





Using integration by parts, we obtain 

E{{vdx,b)) = J {vxfdx + b^. (4) 



In the Euler equation we have A~^m, but our inertia operator A 
is degenerate on vir since A has a non-trivial kernel consisting of con- 
stant vector fields on . It is necessary to consider the quotient space 
Vir /Rot{S^) of the Virasoro group by the subgroup of rotations of the 
circle. This is not difficult since it is sufficient to consider m only from 
the image of the inertia operator, see [Q] for details. 
We obtain the Euler equation on vir* 

^{u{dxf,a) = -ad^-i(„(d^)2,^)(tt(da;)^a) = 

= ~{(2uA^'''Ux + Uxf^^^u — a!v^'^Uxxx){dx)^ ,Q). 
If we put V = Ar^u we obtain two equations: 

ai = 

(so a is a constant) and 

'^txx — ^VxVxx '^'^XXX Ci'^XXX' (5) 

After the change of variables u{x,t) — v{x,t) -I- a we obtain the Hunter- 
Saxton equation 

It has also been shown [Q] that the Hunter-Saxton equation can be con- 
sidered as an Euler equation on the homogeneous space Diff+(S'^)/Rot(5'^) 



5 



of the group of orientation-preserving diffeomorphisms of the circle mod- 
ulo the group of rotations of the circle. The construction is quite analogous 
to the construction on the Virasoro group. It is necessary to consider an 
inertia operator A : vect{S^) — > vect{S^)* given by 

A{vd^) = {-Av){dxf. 

This inertia operator will give an equation 

"^txx — '^Vx'^xx "^yxxxj 

i.e. the Hunter-Saxton equation. This coincides with equation for 
a = 0. We will see that in the discrete case the situation is analogous. 

3 Discrete Lagrangian systems on the 
Virasoro group 

Let M be a manifold, let L be a function on M x M. A discrete Lagrangian 
system describes stationary points of a functional S — S{X) defined on 
the space of sequences X — {xk),Xk G M, fc G Z, by a formal sum 

=^i(a;fe,a;fe+i). 

fcez 

The function L is called the Lagrangian. 

Let us assume that M is a Lie group G. Our goal is to investigate 
potential candidates for a time-discretization of the Euler equation. Since 
the Euler equation describes geodesies on a Lie group equipped with right- 
invariant metric, it is natural to require that the function L is right- 
invariant 

Vg G G L{xg,yg) = L{x,y). 

Let H{x) — L{x,e), where e is the identity element of G. A right- 
symmetric Lagrangian L{x,y) is completely determined by H. Indeed, 
L[x,y) = L{xy-^,e) = H{xy-^). 

We shall consider a class of functions of the type 

Hiif,F))^F^+ Juif'{x))dx, (6) 



where / is a diffeomorphism, F G R. It is a quite natural class since the 
scalar product in vir defining the Hunter-Saxton equation (^) depends 
only on a derivative of a vector field. 
Let us consider a functional 

fcez 

where {{fk,Fk)} is a sequence of points on Vir and L is defined using the 
function H 1^ as described above: 
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Theorem 1 The discrete Euler- Lagrange equations s{/ffF^) ~ ^ '^''^ 
following: 

- Qk + ^k+i = 0, (7) 
+ 2f7.+i (log((..,-^,)'))" + U" (-4rT] = 0, (8) 



v«i)7 («i)')^ 

where (ujk,^k) and {cj^+ij^k+i) ctre discrete analogues of angular veloci 
ties, 

{uju^i) = {fi-^,Fi^^)o{fi,Fir\ I G Z. 
Proof. It sufficient to consider a variation of the form 

{id^ev,eA) o {fk,Fk) 
and to find a Taylor series of tlie variation of 

L((/fc_i,Ffc_i),(/fe,Ffc)) + L((/fc,Ffe),(/,+i,F,+i)) 
up to order O(e^). We obtain 

(id + ev, eA) o (A, 
+L{[id + ev,eA)o{fk,Fk),[fk+i,Fu+^)) = 



-2Q.k (log(ajfe))" + U" {u'k)u'ku'k - 



i-^kli)" 



vdx + 0{e ) =0 



.K+i)7 (K+i)')^ 

Since A is an arbitrary constant and v is an arbitrary periodic function, 
we obtain our formulae for the discrete Euler-Lagrange equation (0,^). □ 
The first equation (Q) is quite simple: it says that Q,k+i = i^k, so Qk 
is an integral of our discrete Lagrangian system. The second equation is 
quite complicated. It is necessary to remark that this is not a differential- 
difference equation. Indeed, it includes tJjT^]^ which is a diffeomorphism 
inverse to Uk+i- So, it is a complicated relation between Uk and Uk+i, and 
it is better to view it as a correspondence (i. e. multivalued mapping) 

LUk + 1 

from Difr+(5^) to Difr+(5'^). 

Let us now find a continuous limit of our equations (|^j^) . Our defini- 
tion of a continuous limit is the following. Firstly, we suppose that the 
angular velocity is of the form 

{LOi,ili) = {id + evi{x),eAi), 

i. e. the angular velocity is the identity element of Vir up to 0(e). Sec- 
ondly, we suppose that 

Vkix) = v(x,t), Ak=A{t), Vk+i{x) = v(x,t + e), Ak+i=A{t + e). 

(9) 
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This is a quite natural discretization of time. We substitute these formulae 
in the Euler-Lagrange equations (QJ^ , find a Taylor series, and finally we 
define a continuous limit as an equation arising in this Taylor series as the 
term of lowest order with respect to e. 

Theorem 2 Let U"Jl) 7^ 0. Then the continuous limit of the Euler- 
Lagrange equations flSj is the following: 

Q — const, (10) 

'^txx — '^Vx'^xx ~1~ riVxxx ^ ^xxx • (-^-^j 

Proof. We substitute the formulae (0) in the equations and we use 
a Taylor series. From the equation (0) we obtain 

thus A{t) is a constant. Let us denote it by Q,. This gives us equation (px|). 
From equation we obtain 

e^{-AQ.Vxxx{x,t) + 2U" (l)vxix,t)vxx{x,t) + U" {l)vVxxx- 

-U"il)vtxx{x,t)) + Ois')=0. 

Since U" (l) 7^ 0, this gives us equation (pT|). □ 

If we do a simple change of time t i—> —t and change our constant 
a = —jp^ij, we obtain the Hunter-Saxton equation exactly in the same 
form as this equation arises as the Euler equation (^): 

'^txx — '^^x'^xx '^^xxx OVxxx- 

As we see, a wide class of discrete Lagrangian systems has the Hunter- 
Saxton equation as a continuous limit. We have a natural question: how 
to find a "correct" discretization? We will discuss it in the next section. 



4 A candidate for a natural integrable 
Lagrangian time-discretization of the 
Hunter-Saxton equation 

Since the Hunter-Saxton equation is integrable, it is natural to suppose 
that a "correct" discretization is integrable in some sense. 

To find a function U giving an integrable discretization we can use 
the observation by A. P. Veselov and J. Moser mentioned in the 

Introduction. Namely that known integrable Lagrangian discretizations 
correspond to right-invariant (or left-invariant) symmetric Lagrangians: 

L{x,y) = L{y,x). 

Thus, let us also assume that our Lagrangian is symmetric. The La- 
grangian is given by the function H : 

L(x,y) = L(2;y~\e) = H{xy~^). 
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It is easy to see that symmetric Lagrangians correspond to inverse- 
invariant function H : 

H(x'^) = H{x). 

Let us consider H of the form (^. Which functions U give us an inverse- 
invariant HI This question was studied in and a sufficient condition 
was found: 



Lemma 1 (l4| // a function U satisfy the condition 



xU 



= U{x), 



(12) 



then a function H : Vir — > R defined by formula (j^ is inverse-invariant: 
H{if,F)-')=H{(f,F)). 
A second derivative of the identity Al2h is 



U" 



U"{x). 



Using this identity we can rewrite equation (n) as 



-2nkiiog({ujky)y -u' 



+ 



+2n,+^{iogiiuj^l,y)y -u' 



(13) 



This is very symmetric expression and it is a complete derivative with 
respect to x. We note that the Hunter-Saxton equation arising from the 
Euler equation (^) is also a complete derivative with respect to x. 

The condition ( p^ is very restrictive, however there are many functions 
satisfying (^) . If we consider other cases studied, there are no indications 
how to find an appropriate U. In the case G = SO{n) g 0, the 
final form of a Lagrangian is obtained by imposing the condition that the 
Lagrangian L{x,y) should be bilinear in x, y. It is not clear what is the 
analog of this condition in our case. In the case G — SDiff(R^) |l2| ^ 
the Lagrangian is chosen as a natural generalization of the SL{2) — Sp{2) 
case. It leads to an integrable system, but there is no evidence that it is 
the unique Lagrangian with this property. 

All we can do is to guess which function U satisfying the condition ( |l2| ) 
gives us an integrable system. The answer to this question is still un- 
known, but we have a good candidate. If we look at the equation (|l^), we 
can see that if we have cuk, we must solve a differential equation to find 
'^fc+i and then use the inversion operation in the Virasoro group to find 
u!k+i- It is natural to suppose that an integrable system corresponds to an 
integrable differential equation for ifZli- least, analogous situation was 
observed in the case of SDiff(R^) [Q. In ||l^ we showed that a function 
U{x) = yfx gives us a system such that if we have lou, we can find a;t;+i 
by solving a linear first-order differential equation and using the inversion 
operation in Vir. In the rest of this section we follow [h 4| . 
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Let us consider the equation ( p^ with U{x) 



1 



-2nfe(iog(K.)'))'- 2^/('^*)'^ 

+2n,+i(log((a.-^,)'))' - yi^kl,)' 
Let us integrate this equation once and put 

^ 1 T 1 



(14) 



We obtain the equation 



, $' *'\ 1 1 
8n -4- + — +- + -+ C = 0, 



where C is a constant of integration and n — Qk = ^k+i (remember that 
fifc is an integral). This equation is equivalent to the equation 

r' r f C 1 $'\ 1 

+ Vjf \ -I =0. 

This is a linear first-order differential equation for "if with periodic co- 
efficients depending on $. For generic $ it has only one solution, so 'I' 
is determined by $ up to a constant C. Reconstructing uj^+i from 'I' we 
obtain another constant, so we have a following result: cuk+i is obtained 
from LJk by a two-parametric correspondence. To find cuk+i starting from 
iOk, we must solve a first-order linear differential equation to find (jJ^+n 
and then reconstruct ujk+i from uj^+i inversion. As mentioned before, 
the analogous situation was observed in the case of SDiff(R2) ||l|] which is 
integrable . For this reason we consider system ( |l^ ) to be a good can- 
didate for an integrable Lagrangian discretization of the Hunter-Saxton 
equation. 



5 Discrete Lagrangian systems on the 
group of orientation-preserving diffeomor- 
phisms of the circle 

As explained in Section ^ we can obtain the Hunter-Saxton equation 
using not only the Virasoro group, but also the group Diff+(S'^), and 
the construction is quite analogous. Similarly, we can consider discrete 
Lagrangian systems on the group Diff+(5''^) and obtain results analogous 
to the results for the Virasoro group. In this Section we formulate these 
results, proofs are omitted since they are analogous to the proofs of the 
last two sections. 

We shall consider a class of functions H : Diff+(S'^) — > R of the type 
H{f)^ J U{f'{x))dx, (15) 
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where / is a diffeomorphism. 
Let us consider a functional 



fcgZ 



where {fk} is a sequence of points on Diff+(S'^) and L is defined using 
the function H ([s]): 



mg: 



Theorem 3 The discrete Euler- Lagrange equation 4j- = «s the follow- 



0, 



(16) 



where ujk and ujk+i are discrete analogues of angular velocities, 

This is exactly equation (^) with fl^ = ^k+i = 0. It is simi- 
lar to the relation between the Euler equations on YiT/Rot{S^) and 
BiS+{S'^)/Rot{S'^), see the end of Section ^. Thus, Lagrangian discrete 
systems on Diff+(5'^) are in this sense the "reductions" of Lagrangian dis- 
crete systems on the Virasoro group. Using this remark, it is easy to find 
a continuous limit of equation ([l^). 

Theorem 4 Let [/"(l) 7^ 0. Then the continuous limit of the Euler- 
Lagrange equation ]ld) is the following: 



Vtx 



Lemma 2 If a function U satisfies the condition then a function 

H : Diff^{S^) — > R defined by formula is inverse-invariant: 

Hif-') = H{f). 

If the function U satisfies the condition (^^, we can rewrite equa- 
tion jl^ ) as 

1 \ . / 1 



U' 



+ U' 



= 0. 



As in the previous Section let us consider U{x) = y/x. Equation (|l 
becomes particularly simple 



(17) 



= 0. 



(18) 



To find u>k+i starting from cuk, we can rewrite the previous equation 

as 

(t.,-^i)' = (C- V^)', 
where C is an arbitrary constant. We see that given Uk, we can find cuk+i 
using integration and the inversion operation in the group Diff+(SM. We 
see that system ( |l^ has the same property as the the equation ( |l^ and 
can be also considered as a good candidate for an integrable Lagrangian 
discretization of the Hunter-Saxton equation. 
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6 Conclusions 



We studied discrete Lagrangian systems on the Virasoro group and the 
group of orientation-preserving difFeomorphisms of the circle. It is shown 
that under some very natural assumptions, these systems have the Hunter- 
Saxton equation as their continuous limit. In particular we studied a 
special discrete Lagrangian system on the Virasoro group which has the 
following properties: 

1. Its continuous limit is the Hunter-Saxton equation. 

2. It can be solved by solving a first-order linear differential equation 
and using the inversion operation in the Virasoro group. 

3. The Lagrangian is symmetric and right-invariant. 

The "reduction" of this system to Diff+(S^) has properties |l| and|^, but ^ 
is replaced but simpler property: it can be solved by integration and the 
inversion operation in the group Diff+(S'M. 

In the case of SDiff(R^) studied in [ p^ , Q the Lagrangian discrete 
system is integrable (in the sense that the dynamics are linearized). We 
conjecture that the same is true for two systems studied in this paper 
since their properties are analogous to the properties of the system on 
SDiff(R2). 

The Lagrangian discrete system on SDiff(R^) can be interpreted as 
a chain of Backlund transformations for the (integrable) Monge-Ampere 
equation ^| . This is nowadays a standard way to think about integrable 
discretizations. It would be very interesting to prove similar theorems for 
the systems discussed in this paper. 
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